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(^ . Abstract: The integrability of the deformed quantum elHptic Calogero-Moser problem introduced by Chalykh, 
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Feigin and Veselov is proven. Explicit recursive formulae for the integrals are found. For integer values of the 
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^ ' parameter this implies the algebraic integrability of the systems. 
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. 1 Introduction 

> 

^O . Deformed quantum Calogero-Moser (CM) systems were introduced by Chalykh, Feigin and Veselov QlEl, who 

O ■ 

\^ ■ proved their integrability in rational and trigonometric cases and conjectured that the same is true in the elliptic 



case. The aim of this paper is to prove this conjecture. 



^p ' Elliptic deformed CM system corresponds to the following Schrodinger operator 

'■^^ " " '■2 



771, 



where all but one "masses" are equal, mi = mT^ , 7712 — ... = ?n„ — 1, 771 is a real parameter, pj — i-^, 
j — 1, . . . , 7i, and p is the classical Weierstrass elliptic function. The case when m is integer is a special one: in that 
IT^ I case a stronger version of integrability (the so-called algebraic integrability) was conjectured ,2,. The first results 

M ' in this direction were found in 131 , where it was proven in the simplest non-trivial case tt, = 3, 777 = 2. 

The main result of the present paper is an explicit recursive formula for the quantum integrals of the elliptic 
deformed CM system. This proves integrability of the system for all 77 and 777 and due to a general recent result by 
Chalykh, Etingof and Oblomkov 4 this also implies the algebraic integrability for integer values of the parameter 777. 

As a by-product we have also new formulae for the integrals of the usual quantum elliptic CM problem, which 
was the subject of many investigations since 1970s (see in particular [S], [Hj, [3IH], pi llOlITT] '). We will be using 
some technical tricks from these important papers. In trigonometric and rational limits we have the formulae for 
the quantum integrals of the corresponding deformed CM systems which are also seem to be new. 



2 Preliminaries and main formulae 

Quantum Hamiltonian of the deformed elliptic CM problem has the following form 

n 

H=-{md^,+dl + ... + dl_,+dl)+2{m + l)J2pik + 2m{m+l) ^ p,fe, (2) 

fc=2 2<j<k<n 

where di = ^- , pjk = p {xj — Xk) ■ Here p is the classical Weierstrass elliptic function 1 2 , which can be determined 
by the differential equation 

(p' iz)f = 4 (p (z) - ei) (p (z) - 62) (p (z) - 63) = 4p3 (z) - 52P (z) ~ 33- (3) 

The Laurent expansion of p at the origin is of the following form |12j 

00 
p{z)=z-'+J2j,kz'\ (4) 

where 

72 = -^52, 74 = -^53. 

The coefficients ^2k are related to the so-called Bernoulli-Hurwitz numbers BH (k) |13j 

1 BH{2k + 2) 
^^'' ~ (2A0! (2fc + 2) ■ 

There is a recursive formula which allows one to obtain the coefficients 72fc+2 from the coefficients of the lower 

order: 

fe-i 

^^^+^^ (fc-l)(2fc + 5) g^^^-^^^-^^' ^ = 2,3,.... (5) 

The relationship lO is easy to verify. One needs to differentiate (PJ to obtain 

p" (z) = 6p2 (z) - 1^2, 

use the expansion Q and collect terms at the appropriate degrees of z. 

To construct the integrals of the operator ^ we will follow the idea going back to jH], [7]. Namely, the 
integrals are constructed from the highest one by successive commutators with some function which in our case is 
m~^xi + X2 + ■ ■ ■ + Xn- This highest integral is of order n and will be denoted below as /. The rest of this section 
is to explain the main ingredients of the formula for /. 

Let us introduce the following differential operators J)*"' in di with constant coefficients 

[5] 
^''^Po,k^'l + Y,P2^,k^t^'. (6) 

i=2 



The constants po.fc , A: = 2, 3, ... , are given by 

Po.fe = ^ nfJi' (1 - M , A: = 2,3,..., (7) 

and constants P2i,ki i — 2,3, . . . , k — 2i,2i + 1, . . . , are determined by the foUowing recursive relationship 

P2i,2i = 0, 1 = 2,3,..., (8) 



i-l 



P2t,k = k_2-i' -P2i,k-1 - {l+m)J2 (2« - 2j)!Cj._2]+i72j-2jP2j-2,fc-l, 

i = 2,3,..., A: = 2i + l,2i + 2, ... 



(9) 



and C^ = ^|.^_^^s, is a binomial coefficient. To illustrate the formula © let us write explicitly the formulae for the 
first few values of i : 

_ (l-m(fc-5)) ^ /-, , \ nLi(fc~0 „ „ l, _ c, G 

P4,fc — iP4 P4,k-1 — (i^ + m) 3 72P0,fc-l, K — 0,D, . . . , 

pe.k = <^'T!V'" P6.fc-i - (1 + m) "'-f'S 4Po.fc-i, fc = 7, 8, . . . , 

P8,fc = ^^'T!8~'» P8,fc-i - (1 + m) ( n'=f-S 6Po.fe-i - %72P4,fc-i) , fc = 9, 10, . . . . 

It is interesting to note that for the special values of parameter m = j, where I is a positive integer number, most 
of the constants Q and are zero. For example, if to = i only po 2 is non zero, if ?7i = i then only po,2 and po 3 
are non zero, if to = i then only po,2, P0.3, Po,4 and ^4,5, p^^e = Po,2P4,5, P4,7 = Po,3P4.5, P4.8 = Po,4P4,5 are non zero 
and so on. 

Let us introduce the following notations: 

<;j = (to + 1) <r {xi — Xj) , 2 < j < n, 

uij = (to + 1) p {xi — Xj) , 2 < j < n, 

Uki ^ m {m + 1) p {xk " xi) , 2 < k <l <n, 

where <i is the standard elliptic (^-function: ''j^' — —p{z) . 

We will need also to consider all the subsystems of the deformed CM system. Let S" be a subset of the set 
{l,2,...,n} and a = {ii, J2, • ■ ■ , jt} , ji < J2 < ■ ■ ■ < it, be a subset of t indices chosen from S. The set of all 
different subsets a of size t of the set S will be denoted by 

6 (S*; t):::^ {a = {ji,J2, ■■■Jt} ■■ ji < 32 < ■■■ < jt, jk e S} . 

If (T G S {S; t) define the set S\(j = {j : j <E S and j ^ cr} . If S* = {1, 2, . . . , n} we will use short notation & (t) — 
& ({1, 2, . . . , n} ; i) and ci = {1, 2, . . . , n} \cr. If a contains only one element a = {k} the brackets will be omitted: 
k will denote a set which contains one element {k} , S\k — S\ {k} and k = {1, 2, . . . , 71} \ {k} . We will also use 
notation aia2 to denote the intersection of the subsets cti and (72 : (TiO'2 = (Ti n (T2. 



We use the notation ad*^ (S*^) to denote the repeated commutator 

a<(D^) = [q„...[q,,[q,,S)'^']]], a^{ji,j2,...,jt}. (10) 

Note that the order in which q^. are used is not important because of the form of the operator J)*^. Define 

*-i(TGS(i;t) *-0(Tee(i:t) 

Let the set 5* consist of k elements and contain 1, then define the operator Qg as 

[I] [I] 

0^ = ®'+E E <(^'"*)=E E -'^K^^'^^O- (11) 

t=lcree(S\l;t) t=0 cre6{S\l;t) 

We also use the notation ds to denote the product ds = Yi ^^j- ^^ ^ ^ ^ then A5 denotes the sum A5 = 

"^^1 + X^jesXi ^|- If 1 ^ "S* then As = X^jes ^|- ^^^ ^-S ^^ '^^^^ mean the corresponding quantum integral of the 
system with the Hamiltonian 



Hs^'As + 2j2 



Ujk- 



j<k 
i,kes 



Now we are ready to give the formula for the integrals; 

n-2 



^=E(-1)*^' E hd, + {-ir{n-l)di...dn+X, (12) 

t=i cree{t) 



where 



^ = ®+ E E ^-Q- (13) 

*=1 cre6(i;2t) 

and Xa- are related to a non deformed CM subsystem and are determined by the recurrent formulae 

Xi = J2'j=2 ^{j:n}^i\{j,n}^ if n = 2p - 1; Xj = 0, if n = 2p and X{2,3} = "23- (14) 

Theorem 1. The operator I defined by |Ji^) and <J5^) commutes with the deformed elliptic CM operator H . 
Remark. The formula (|12|) is valid in the non deformed case also: the operator 

■;i-3 

A = E(-1)*^' E Iad, + {-^T{n-2)d2...dn+X-, 
*=i (Tee(i,t) 



commutes with the operator 



2<j<*:<n 



which is the usual n — 1 particle elliptic CM operator. 



Idea of the proof of Theorem 1. The proof will be done by induction. The main idea behind formula (|12|l 
consists in the observation that one can use the commutativity of lad^ with H^ to simplify the commutator [/, H\ 
to the expression 

[/, H] = [X, H] + Y^'' \x,d,,2 Yj=, u,} + Ei<,<,<„ [X^d,di,2u,i] , 

where terms X, Xj and Xj^i depend only on di and Usr, 1 < s < r < n. This is shown in Lemma 1 in section 6. At 
the next step we notice that if X is given by H13() the commutator [/, H] can be simplified further to the following 
expression 



i=2 



0i5,,2^ 



U-jl 



1=1 



^ [Qj,i;,2{uik + Uu)]X{kd} 



2<k<l<n 



where Q, Qj and Oj.^- depend only on di and Mis, 2 < s < n. From this we can deduce that 

^® [ n 1 

and, therefore, it seems natural to use the operators ad* to construct 8. At this stage the only freedom left is 
in choosing operators J)'"' which must be the operators in di with constant coefficients. To ensure that [/, H] = 
these operators must satisfy the relation 

[D",p(xi-a;,)] + (l + TO) [[<r(xi-2:,),S""^] ,p{xi-x^)] 

+m [D"-i, p {x, ~ X,)] d, - ^-^p' (xi - X,) S)"-i - (i^2)"-ip' (^i - x,) = 0. 
which is equivalent to a large set of identities. It is remarkable that the constants in J)*^ can be chosen in such 
a way that all these identities is satisfied (Lemma 3 of Section 6). The choice of the constants is related to the 
BernouUi-Hurwitz numbers and is described above. The complete proof of the theorem is quite technical and is 
given in a separate section. 

3 Examples: formulae for two, three and four particles. 

To illuminate our formulae let us consider more explicitly the case of small n. 



3.1 Two-particle case. 

In that case we have the operator 



H = -mdl -dl + 2{m + l) pxi = -mdl - d^ + 2ui2 

which is trivially integrable since the operator di + 82 obviously commutes with it. This system gives the formula 
for operator /{1.2} which starts the recursive construction of the integral / (|12|l . We have 

/ = i{i.2} = ^[H + idi + d^f) = did2 + ii^a2 + {m + 1) P12 = 8182 + 2)2 + wi2 = 8182 + e, 

^ = ^{1,2} = ©{1,2} =D2 + ui2 = S)2 + [<r2,2)]. 



3.2 Three-particle case. 

The integrals for the problem of three particles have been found in 14 . The operator of the third order has the 
form 

+ ii:_!Iil(m + 1) ((pi2 + P13) 5i + 9i (pi2 + P13)) 

= 9192^3 + P0,2(92 + 83)01 + U23dl + Uis^a + 1*1293 +&^+ [c^2, 3^] + [<r3, 2)2] , 

Operator / can be rewritten as 

I = I{2..3}dl + I{l,3}d2 + I{l,2}d3 - 2818283 + X, 

where 

X = e = J)^ + [<r2, S^] + [<^3, 2)^] , and /{2.3} = ^293 + U23. 

3.3 Four-particle case. 

One can check by direct calculation that the operator 

/ = 9i929394 + ^i^ (9293 + 9294 + 9394) 9f 

^ (l-mKl^2m) ^Q^ ^Q^^ Q^-^ Q3 ^ (l-r»)(l-2m)(l-3m) g4 

+TO(m + 1) (p349i92 + p249i93 + p239i94) + (m + 1) (pi49293 + P139294 + P129394) 

^ (1-m) ^^ ^ -^^ Q^ j.j.^^_^ ^ ^^^^ Q^ ^ j-^^^ ^ ^^^^ ^^ ^ j.^^_^ ^ ^^^^1 ^^^1 
^ (1-m) j.^^ ^ -^^ j.|-^^^ _^ ^^^^ ^_^ ^ ^^^^ ^ ^^^^ ^^ ^ ^^^_^ _^ ^^^^ g^-j Q^ 
^ (l-mKl-2m) j.^ _|_ -^^ ^^^^_^ ^ ^^^ _^ ^^^^ g2 ^ g^ (^^^ _,_ ^^^ _,_ ^^^^ g^ _^ q2 (^^^ _^ ^^^ _^ ^^^-jj 

+ (i_!Ii)rri(m + 1) (P34 + P24 + P23) 9? 

+ (1 - m) (m + 1) (pi4p24 + Pi4p34 + P24P34) + m (m + 1) (pi2p34 + Pi3p24 + P14P23) 

commutes with H. The operator / can be rewritten in the following recursive form 

/ = 39i929394 + X + /{2,3,4}9l +-^{l,3.4}92 +^{l,2,4}53 + ^{l,2,3}i94 (15) 

~-^{3,4}l9lt^2 — -f{2,4}f^ll93 — ^{2,3}^lt'4 — ^{l,4}t?293 — /{i,3}9294 — /{l,2}t^394, 



where 



X = 6 + X{2,3}0{1,4} +-'^{2,4}0{l,3} +-''^{3,4}0{l,2} 

= 1)4+ [^2,2)3] + [^3,1)3] + [.^4,1)3] + [^2, ['^S,®']] + ['^2, [q,®']] + [<?3 , [^4,2)2]] 
+U23 (S^ + ^4, S]) + U24 (S^ + ^3, 2)]) + U34 (S^ + ^2, ©]) • 



and 

-^{2,3,4} = did2d-i + U2i,di + ui^d2 + W1293. 

4 Integrability of the deformed elliptic quantum CM problem. 

Let us introduce the function 6 = m^^xi + X2 + ■ ■ ■ + Xn and consider the corresponding "ad" -operation: 

adg{L)^[9,L]. 

Following to the procedure known for the usual CM system (see, for example, |S]) consider the operators Lk = 
adg (/) , A: = 0, 1, . . . ,n - 1. 

Theorem 2. Operators Lk = adg (/) , fc = 0, 1, . . . , n — 1, where I are given by 1^1S9\) commute with each other 
and with operator H. 

Proof. The proof is similar to the non-deformed case |S1 E] . 

Let us first prove that Lk commute with H and with X]7=i ^i- Proof is by induction in fc. For fc == it follows 
from the Theorem 1. Let us assume that this is true for k — i, then for k ~ i + 1 we have by the Jacobi identity 

and 

= [{m^^+n-l)ld,L';']+[0,e]^0. 

To prove that operators Lk and L;, k y^ I, commute one can use the arguments of the paper |11| . Consider an 
involution S on the space of all differential operators on i?" corresponding to the change x -^ —x and the standard 
anti-involution *: operator L* is a formal adjoint to L. We have [L'J,L2] = [-^11^2] and [L5;,i2] — ~ [iij-^2]* ■ 
Operators Lk have the following properties with respect to these involutions: L^ = Lf = (—1) Lk- Now, consider 
the commutator C — [Lk, Li]. By the Jacobi identity [C, H] ~ 0, therefore we can use Berezin's lemma ^S] which 
states that in such a case the highest symbol of C is polynomial in x. In this case it is also periodic, hence the 
highest symbol must be constant. Now, we also have 

C' = [Lk,Lif = [LlLf] = [LIL:] = - [Lk,Lir = -C* . (16) 

Since the highest symbol of C is constant the highest symbols of C* and C^ are the same but it follows from (|16(l 
they they must be different by a sign. This means that the highest symbol of C is zero and hence C is zero. 

Theorem 3. Deformed quantum CM problem 0) is integrahle for all n and m and algebraically integrable for 
integer m. 

Proof. The complete family of the commuting quantum integrals for arbitrary m is given by the previous 
theorems. The algebraic integrability in the case when m is integer follows from the general result due to Chalykh, 
Etingof and Oblomkov (see Theorem 3.8 in ^). 



5 Trigonometric and rational degenerations 

Trigonometric degenerations of the Weierstrass p-function corresponds to the case when one of the half periods 
LUi or UJ2 is infinite, which happens when two of the roots of the polynomial (Q collide. For example, the case of 
ei —62—0, and e-? — —2a corresponds to uj — cxd, uj = i—^= and p (z) — a A — . , .?°, Choosing a = ^ we have 

Vl2a c \ / sinh^ vSaz ^ -^ 

f'i^'' - 3 + ;hIT?7 -^ ~Z^fc=i l2iq:2j-(2l)r^ ' <rizj - -3^ + cotti2;, and 72k - -j2k+^-(2i~)r^ 
where B2k+2 are the classical Bernoulli numbers defined by the expansion 

z — 1 _ 1 7 4_ V°° ■B2fc 72fc 
e'-l ~ 2'^ '^ 2^k=l (2ky. ' 

In this case the Hamiltonian H takes the form 
H = - {md! +di + ... + dl_, + dl) + (=±iKli^ (1 - rn + ^f) + Y^^^^ ^i^^^, + E.<,<.<« s^^^^Ly 

Therefore, the formulae for the integrals in this case can be obtained using the following recursive formulae for 
constants P2i,k ■ 

P2i,k = '- )rj2l -P2^,k-l + (1 +"l) L '-^fc-2j + l (2»-2j+2) -B2»-2j+2P2j-2,fc-l ■ 

j¥2 

The integrability of this problem was shown in 2 . In 11 6| a recurrent formula was found for the quantum integrals 
with the highest symbols given by the deformed Newton sums. Our formulae correspond to the deformed elementary 
symmetric polynomials and seem to be new even in that degenerate case. 

The rational degeneration corresponds to both periods be equal to infinity. In this case p (z) = z^'^ and all 
l2k = 0- Therefore, in this case only constants po.fci k — 1,2,..., are non zero. 

6 Proof of Theorem 1. 

We prove Theorem 1 by induction in n. For small n we showed this in the section 3. Now assume that the statement 
of the theorem is true for all fc < n and show that it is true for k — n. 

Under this assumption, let us first show that commutator [/, H] can be reduced to an expression on the 
additional terms X, Xj and X^^. 

Lemma 1. 

[I,H] = [X,H]+jy^^^ \xsd„2Y,l,u,] +Ei<,<;<„ [XM9kdi,2u,i] . (17) 

Proof. Using p2|) we obtain 

Ti-2 

lI,H] = Y,{-lf+' Y. [hd.:H] + i-^r{n-l)[di...dn,H] + [X,H], 
t=i <jGe(t) 



which is equal to 



E 



l3d,,Hs-A,+2Y, 



Ujk 



fe=i 



E(-ir^ E 

t=2 (Tee(t) 



/^a,,i/^-A, + 2^^w,fe-2 Y, 



Ujk 



j£a fc=l 



j<k 

j,ke<T 



+ (-ir(n-l) 



l<j</£<n 



[X,i7] 



Since /^ commutes with H^ (by the induction assumption) and with — Ao- (since it does not depend on Xj, i G cr) 
we have 






-l 




- 


- 


n 

E"jfe 
fc=i 


t^2 aeeit) 


n 

/^9^,2^^-Ujfe - 

jea fe=l 


- 2 ) ; Ujk 

j<k 

],kea 



+ (-ir(n-l) 
Now we use (|12() again and obtain 



l<j<k<n 



[X,H]^ 



ii^H] = y: 

3 = i 



n— 3 






E 



(-1)""' (n - 2) ai . . . 9„ + Xj9,, 2 ^ Ujk 

fe=i 



n-2 



E(-ir^ E 

o-Ge(t) 



t=2 



^*5-25]E"^'^-2E 



Ujfc 



j^(T k—1 j<k 



+ {-lT{n-l) 
If we cancel the repeated terms we obtain 



di...dn,2 ^ Ujk 

l<j<k<n 



[X,H]. 






{-ir-'in-2)di...dn+Xjdj,2j2 



Ujk 



fc=i 



n-2 



E(-i)' E 

i=2 (je6(t) 



Iad„,2 ^ MjTi 



]<k 

j,k£a 



+ {-lT{n-l) 



di...d„,2 Y ^jk 

l<j<k<n 



[X,H] 



We use (|T^ again 






(-1)""' (n - 2) 9i . . . 9„ + X,'d,,2 Y, ujk 



fe=i 



E 

l<j<k<n 



*=i CTee(ifc;t) 



n-2 



i<i<fe<" *=3 cree(t) 



■ 




/*9, 


2E 



Ujk 



+ (-ir(n-l) 
Finally, canceling the repeated terms, we get 



l<j<k<n 



[X,H]. 



[I,H] = (-ir2 



Si .. .a„,- (n- 2)^^Ujfc + (n-3) ^ u^/c + (?i - 1) ^ Ujfe 

l<j<fc<n l<j<fe<n 



j=l fc=l 
Mi 



J=2 



fc=l 



E [ 



Xff.djdk,2ujk 



l<j<k<n 



which simplifies to 



[I,H] = [X,H]+Y 
i=i 

Lemma 1 is proven. 

Lemma 2. In the non deformed case 

n 
J=2 

Proof. To prove that 



^i5j>2^Mji 



fc=i 



+ E [%5jafc,2Mj7c 

l<j<fe<n 



Xya„2^ 



Wj; 



i=2 



E [^i«A.a,,2ufe,]=o. 



2<fc<i<ra 



i=2 



Xy9„2^ 



W,7 



i=2 



E [Xuidkdi,2uki\ 



2<k<l<n 



(18) 



one must repeat the arguments of Lemma 1. We will need the addition theorem for the Weierstrass elliptic 
function ^21 

/ Pij Pjk pki \ 



Tijk = det 



Pij Pjk Pki 

V 1 1 1 y 



= 0. 



10 



Two cases must be considered to prove the lemma: n is odd and n is even. If n is even (JTS)) is reduced to 



n n r 



2EE 

j=2 1=2 

1^3 



Xj^.dj,Uji 



0, as Xj = and X^j^^ = 0, 2 < k < I < n. In this case we have 



n n 



j=2 1=2 
1^3 



■11 7i /(. tl 

= T.T. ^ii [^^- ,",/]= ™ (m + 1) ^ ^ Xi -p^, 

i=2 i=2 j=2 1=2 

1^3 1^3 

. (m + 1) Y. ^3i (^ii - ^if) = "^ (™ + 1) E ^^J'^^ 



= ?Ti 



iijk 



2<j<l<n 



2<i<j<k<n 



If n is odd lfTH|l becomes [Xj, — Aj] + ^ [Xj^j-3fe9;, 2ufe;] = 0, since in this case Xj- = 0, j = 2, . . . , n. We 

2<i:<i<n 

have 



[X 



2<k<l<n k=2 ^ k=2 '^ 2<k<l<n 



9k - uli) 



E^-2 E ^i^.-<. + 2E i|-E^u.-H'^ = '- 

fc=2 ^^fe 2<fc<i<n fc=2 I '^^'^ 1=2 j 

~ ~ \ l^k / 



= E 

fc=2 

Lemma 2 is proven. 

Lemma 3. The operators J)" satisfy the relation 

[D",p(xi-a;0] + (l + m) [[<; {xi - x^) ,D''-^] ,p{xl~x^)] 

+m [D"-i, p (x, - X,)] a - ii±^p' (xi - X,) D"-i - ii^D"-ip' (xi - x,) = 



(19) 



for any i — 2, . . . ,n, where d — ■^—. 
Proof. 

Let us denote the left hand side of the relation (|19|l by 3^„ . 
It can be shown by a simple direct calculation that H19() is true if n = 1, 2, 3, 4. The equations in these cases are 

n = 1 : 3^1 = [9, p] + m [1, p] 9 - (i±^p' - ii^p' = 0, 

n = 2 : 3^2 = [i^52, p] + (1 + m) [[,, 9] , p] + m [5, p] a - iii^p'a - ii^ap' - 0, 

r l-2m . 



id\p\ + (1 + m) [[^,92] ,p] +rn [92,p] 9- ii±=^p'92 - ii^92p' = 0, 



„ = 4 : 3^4 = [i^94, p] + (1 + m) [[.;, 93] , p] + m [93, p] 9 - ii±=^p'93 - ii^9V = 0. 



For n > 5 we use © to obtain 

yn = [po,„9", p] + (1 + to) [[^,po,n-i9"-i] ^ p] + ™ [po,„-i9"-\ p] 9 - ^^^p'po,„_i9"-i 



-ii^Po.«-i9"-V' 



Ep2,,„9"-2\p 

i=2 



+ (1 + to) 



+771 



Y. P2i,n-ld 

4=2 



in — 1 — 2i 



(1+m) , / 



P' Z) P2i,«-l9 
i=2 



<;, J2 P2i,n~ld'' 
i=2 



in-l-2i (1-™) Y^ „ C)n-l-2j,„/ 

4=2 



'P' 
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We can express po,n through po.n-i and P2i,n through P2i,n-i, I — 0, ■ ■ ■2i, a.s described by lO and obtain 



yn = 



'~"l"''^ Po,n-ia", pi + (1 + m) [[^,po,„-i5"-i] , p] + ^ [po,„-i9"-i, p] Q 



(1+™) .„'^ r)"-l (!-'") ^ an-1,^/ nj-Tn"; 
2 PP0,n-lO 2 PO,n-lC p — (_i+mj 



1=2 



Z-/ ,i-2j Wi,n-lC^ , P 

2 



+ (1 + to) 



■^7 Z) P2i,7l^ld 
i=2 



\n—l — 2i 



%P 



O 2 P ^ P2i,Ti-lC 2 1^ P2i,n~lCI p 



-(1 + m) E 

i=2 



E (2j - 2* - 2)!C;:Z2'/_i72,_2,-2P2.,„-i9"-2^ p 

j=t+2 



Denote 



Wk = 



then 



and for fc > 5 



i-m(/c-i) gfc^ p] + (1 + m) [[<;, a'-^-i] , p] + m [d'^-\p] d 



3^7 = P0,6 Wt + P4,63^3 + Pe.e^l , 3^8 = PO.T^S + P4,73^4 + P6,73^2 



fc-3 



3^2fe-l = P0,2fc-2W2fe-l + 2^P2i,2fc-2VV2fe-l-2i + P2fc-4,2fc-2y3 + P2fe-2,2fe-23^1 , 

i=2 
fe-3 

3^2fc = P0,2fc-lW2fc + 2^ P2i,2fe-lW2fc-2i +P2fc-4,2fc-iy4 +P2A;-2, 2/0-1^2- 



i=2 



Below we show that W„ = for n > 5. Firstly, 

^Ws = 2[94,p]9-3aV + 10[[^,9^],p]-5p'a4-80[729\p]+29[93,p]9 

-3p'a4 - 394p' + loa [[,, d^] , p] + 10 [a, p] [,, d^] + lo [c^, d] [d\ p] - so [-/2d\ p] . 

We calculate that d[[<i,d^] ,p] = -ld[d^,p] d + ^dp'O^ + jd'^p', therefore 



10 
l+m 



w. 



-i9 [93, p]d- ip'd^ - \d^p' + hdp'd^ + lOp' [<;, 93] + lOp [53, p] - 8O72P' 

-5p"'92-|p(4)9_lp(5)_8072p' 

-lOp' (3pa2 + 3p'9 + p") + lOp (3p'92 + 3p"a + p'") = 0. 
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In this calculation the identities which are the derivatives of the differential equation of the Weierstrass elliptic 
function Q have been used. At the next step we use the induction assumption that for any s < k Ws — 0. Then 
Wfc can be simplified into 



m- 



1+ra [fik—l 



[9^-1, p]d- ii±!ii)^a''-V' + (1 + m) [[,, a'^-1] , p] 



.(l+m)^,gfc-l 



J=2 L 



We consider 



lik^2l 



jl^Wk = i[9'=-i,p]5-*^9'=-V'+[[<^,a''-i],p]-ip'9^-i 



j=2 L 



T^W,, = i [a^'-\ p] 9 - ^9^'- V + [b, d^^-^] , p] - ip'92'-i 



i-i 



J=2 



E (2j--2)!C2^;:f^ 72,-23 



2l~2] 



,p 



and if fc = 2Z + 1 



I+;;7>V2/+i 



(5ITT) [9^ P] 5 - 4md''P' + [^' ^^'] ' P] - Wd' 

j=2 L 



We show below that yV2/ = 0. The case of W21+1 can be dealt with in the same way. 



(-1 r 



i=2 



ip'92'-i „ E (2j - 2)\Cl\-_l'i2j^2d^'^'^. 



P 



= ^ [992^-2, p] a _ 1^9a2/-2p, + [[^^ 9g2/-2] ^ ^] 



I-l 



-ip'^a^'-^ - E (2j - 2)!C2',':i^^72,-2aa2'-2^-i, p 

J=2 L 

^9 [a2'-2, p] a _ i^dd-^'-^p' - i^p'a92'-2 + a [[,, a2/-2] 



l-l r 



-a E (2j - 2)!c|;:f 72,-2a2'-2^-\ ^ „ ^ (2^. „ 2)\cl-_l' ^12,-29' d^^-^^-\ 

J=2 L J J=2 



/-I 



From the inductive assumption we have that W2i-i = 0, then 

^^'-^^PJ = - 



[^'^''"']'p] = -2I^[a^'-^p]5+2fE^a2'-v+^P'9^ 



-E (2j-2)!c,^;:^2^72,-2a2'-i-^^P 

i=2 L 
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and therefore 



-1-^W2/ 






i=2 
i-1 



(2'-2)! , «2;-l-2j- 

(2;-2j)!72j-2C 



P 



p'[,,a2'-2]-p[a2'-2,p] 



+ E(2j-2)!C|;:f72,-2P'9^'-^^- 

J=2 



2J(2Z 



i^iya [a^'-i, p] - ip"a2'-2 _ ip'a^'-i - p' [,, a^'-^] _ ^ ^^21-2^ ^] 



;-i 



+5 E 7Kt72,-2 [a^'-^-^^ P] + E (2j - 2)!c|;i,^^72,_2P'a^'-^^-^ 



i=2 



(2i-2i) 



J=2 



The commutators can be rewritten as 

2(-2 



21-3 



[d^^'\ p] = J2 ci\-_i-^p^^'-'-'^d\ [9^'-^ p] = y: ci\zi-^p^'^-^-^^d\ 



fc=0 
2(-3 



*;=0 



2;-2-2j" 



[>?,92'-2] = ^ ^2i-2-fc^(2/-3-fc)gfe^ ^Q2l-l-2j^^-^^ ^ ^2i-l-2j-fe ^(2;-l-2,-fe) 5/= 



fe=0 



fc=0 



and, hence, 



2/-1 



2J-2 



T+;;r>^2/ - 2i(2;-i) 2^ ^2/-iP c* + 2i{2i-i) 2^ '-'2i-i P c* 



fc=l 



fe=0 



2;-3 



2i-3 



lp"a2i-2 _ ^p'a2'-l _ ^ ^2J-2-fe^,^(2;-3-fc)5fc - ^ C^'t' " Vp^''"'"'^ 9^= 



2 
J-1 



fe=0 fe=0 

Y^ jzt-zjj^ -f^'" f^2l-2j-k I2l-2j-k)fik 

E (2;-2j)!72j-2 2^_ L.2,_i_2jP ^^ 



2/-1-2J" 



(2i-2)! 
,=2 ^^'-^^'^ 

'' E^'' cl^.^^f^pf^'-^.-.j^fe + 'j2 (2j - 2)!c,^;:f 72,-2P'a2'-^^-^ 



Y^ (2;-2)! 
' 2^ (2;-2j)!T2j-2 
j=2 k=0 



l-l 

E 

i=2 



which can be rewritten as 



-Th^2i 



21-1 



21-2 



1 . V r^l-k,j2l-k)f)k I 1 V (^2;-!-*: (2;-fe)Q, 

T^rry 2^ <^2;-iP " + 21(21-1) 2^ '-"2;-i P c* 



2;(2/-i) 



fe=l 



*;=0 



2/-3 



lp"a2'-2 _ ip'92i-l _ Y: Cllzt'' (P'P^''"'"''^ + Pp(2'-2-fc)) d'' 



k=0 



21-5 



l-l- 



, V- V- (2i-2)! (^2!-2j-fc (2;-2j-fe)a/c 

^^ Z^ Z^ (2;-27)! '2j-2L^2;-l-2j P " 



fc=l j=2 



21-6 

-E 

fc=0 i=2 



Y^ Y^ (2(-2)! 

2^ 2^ (2i-2j)! 72^-2 



1-2 



C^2i-l-T'p^^'-''-''^d''+E i2l-2k-2)\Clti-f2i-2k-2P'd 

k=i 



/Q2fe-1 
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It is shown below that the coefficients by all i9'^, k — 0,1, . . .21 — 1, are zero. We have 

Q2/-2 . 1 /^2 » I 1 f^l ^" _ 1 ^" _ n 

*^ • 2K2;-l)^2;-lF + 2;(2/-l)^2Z-lP 2*^ ~ "' 

2/_3 . 2K2H)C'2Vip(3) + .^C-^ViP^S) _ 2C2V2p'p 

(2;-l)^2/-l + 2;(2;-l)^2Z-l ^ 2C2;_2 ) P P = 
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a2/-4 



K2H) G2ViP<"^ + 2K2H) c-^VipW - ci,_, {p'p' + pp('^) 



When k = 2l-2q-l,q = 2, ...,l - I, coefficient by a2/-2q-i jg 



q-l 



r. , _ (2i-2)! p^^-"^^' _ p'py-'i-->+ppy-'i-^> V. 72,-2 ,^(2q+l-2n , _2q_^ , 

^2l-2q-l — (2;-2g-l)! I (2q+l)! (2q-l)! ^ f^ (2q+l-2j)<.P ^ 2q-l '2(?-2P I • 

To show that this coefficient is zero let us calculate its Laurent expansion and show that it consists of the terms 
by the positive degrees of z only. Since it is also a doubly periodic function it can only be zero. The Laurent 
expansion for the derivatives of the Weierstrass function are given by 

p' = -2z-^ + 0{z), 
pf^-J+i) = -(2q + 2)!z-(29+3)+o(2)^ 

^(2,+l-2j) ^ -(2g + 2-2j)!z-(25-2j+3)+0(;2). 

Therefore, we can find that 

pp^'"-'^ = -(2g)!z-(2«+3)-(2g)!|]72.^-(''+^-'^'+0(z), 

q-l 

p'p(29-2) ^ _2 (2q - l)!z-(29+3) -2{2q- 2)\-f2q-2Z-^ + {2q - iy.J2 2i72»2-(29+i-2») + o {z) 



and 



Therefore 



(.V<---W----') ^ (2, + 2) z-(^«+3) + '^^ (2, _ 2.) 72.z-(^^+^-^^) + 21^72,-2.-^ + O (z) , 

2—1 



^2^-2,-1 = (2!!^-!)! - (2? + 2) ^-(^"+3) + (2g + 2) z-(2«+3) + ^^'f (2g - 2z) 72.^-(^«+^-2*) 

+ 21^72,-2^-^ - Y.T=l 72.-2 (2g + 2 - 2j) z-(2.-2,+3) _ ^i,^_,^^z-^ + O (,) 
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When k = 2l-2q;q^3, ..., I - 1, coefficient by d^'-^i is 

^2l-2q - (^2l~2qy. \llq)\P (2g-2)! + 2^j=2 (2q-^2j)! '2j-2S3 I 

The Laurent expansion for the derivatives of the Weierstrass function are given by 

p(2«) ^ (2g+l)!z-(2^+2) + (2g)!72, + 0(z^), 
p(2,-2j) ^ (2g-2j + l)!z-(29-2j"+2) + (29-2j)!72,-2,+0(z2) 

and we calculate that 

pp(29-2) ^ (2g - l)!z-(29+2) + (2^ _ 2)12(772,-22-2 + (2g-l)!(g+ 1)729 
+ (2,-l)!('x:72.2-(2^-2^')+0(z2), 

pV-3) = 2(2g-2)!z-(2«+2)-(2g-2)!2g72,-2Z-2-(2g-2)!(2g)(^^ 



and 



- ^"''"72.-'2y^'^ = - (2? + 1) --''"^'^ E (2? - 1 - 2^) 72.z-(2^-2^) _ fil±m^\ ^^^ + o (,2 



Therefore 



^2/-2, = fi^ (2g + 1) z-(2«+2) + 72, ~ {2q + 1) z"(2«+2) - E?;' (2g - 1 - 2.) J2^z-^'^-'^^ 



(^i±^f^) 72, + E|;2 72,-2 (2g - 2j + 1) z-(2'^-2.+2) + 1^ 72^_272,_2, + O 



9-1 

E' 

J=2 

«-l 



= (i^ (^-^^^±^^72,+ E 72,-272,-2.+0(z2)j =0. 

Here we have used the identity |SJ- 

The last case to consider is when fc = 0. In this case we need to show that 
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is zero. We have 



Ko = {21- 2)1 \ (2; + l)^-(2i+2)^^2,-(2/ + l)z-(2;+2)_'^'(2Z-l-2z)72,z-(2;-2^) 



(^ 



±^^^) 72. + E 72,-2 {21 - 2j + 1) z-(2'-2.+2) ^ g ^,^._2^,,_2, + O {Z^) 

' i=2 i=2 



{21 - 2)\ (-M(^^, + ^-1 72,-272.-2, + O (.2)) = 0. 



The proof of Lemma 3 is now finished. 

Lemma 4. The following identity holds 



J=2 



,"=2 ,"=2 j=2 k=2 1=2 



l^k 



Here again d = g|-. 

Proof. From Lemma 3 



[S","i,] + [fo,S)""'] ,wi,] +m [l)"^i,uij] a- ii±2i)y'^^.j) 



(!+'") „y Tin-l (l-™)>T-iTt-l 



£(«-lu'^^. = 



Note also that 



E E aC([S"-*-\«i,]9) ^ E E adl{[j)-'*-\u,,])d 
J=2<Te6(ii;t) J=2CTee(ij;t) 



+ E E E.ee(ifef;t) K. (»""''*) '^i'^-] "l^ 

fc=2 i=2 ^ ' 

We sum (|20|) and then use pi|l to obtain 8 and 8^'. We have 



(20) 






(l+"^) „y ^n-i-l (l-m) ^„-t-l„y 
2 "ij"^ 2 "^ "l, 



EE E [a<(S)"-*),«i,]+ E E E 



<{w(®""*"')'"l^- 



1^1 

2 J 71 n 



+™ E E E a<([S)"-*-i,uy])9 + m E EE E K^(2)"-' ^j,iiu|uii 



-2-n 



t=o fc=2;=2CTge(^jfc;-.j) 



■^ E E E <^-aC (D"-*-i) - (i^ E E E aC (»""*"')<, 



t=0 j=2 ^^e{ij;t) 
L 2 J n 

E E E 

t=0 ,=2„ge(y.t) 



t=0 ,=2^gg(y.t) 
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*=o^ee(i;t) ^ 



mil 

2 J n 



"^ E E E adl{[^-'''\u,,])d 



2 J n n 



-^ E E E E [adl{'D"-'-*),u,k]uu 



(l+m) 



' r,At i^n-t-l\ _ (1-m) 



E E E <,<. (D"-*-i) - 

t=o j=2 „^e{i3;t) 



^ E E E K-„(s 



■*^(D"-*-i)u'y 



e, J2 "ij 

J=2 



m J2 [©i, uij] d - ■^^ E wijQj - ^^ E ^jKj + "i E E [©m. wij "n- 



J=2 



J=2 



J=2 



A:=2 i=2 



Lemma 4 is proven. 

Lemma 5. The following identity holds 



[A,e] = 2^[ziu,e^]afe-2mK;][ui,,e^] ai-(i + m)^K„e^]+m^^[z.ife,[ziu,e^,-]]. (21) 



fe=2 



Vfe=2 



fe=2 



fe=2 /=2 



Proof. 



[A,e] = 2™— a. + 25:— a. + n.^ + g^. 



9xi 



fc=2 



fc=2 



^^ 2 2 

From ini) it follows that §^ = [u.k.Qk] and |^ = - ^ ["ife,0fe] , therefore |;§ = - Kfe,eJ and fj 

; — o fc 1 



/c=2 



^ E ["ife'Qfe] + E E b^ik, [""!/, 0fc[]] ,and we arrive at (jnj- 

fe=2 fc=2 ;=2 

Lemma 6. 



J=2 



^jdj.^Y, 



Ujl 



1=1 



J2 [e^.[,2(uu + MiO]X{fc^;}=0. 



2<k<l<n 



Proof. We have 



R=[A,Q] + 2J2 [0, nik] + 2 X^ E [%' "J'] ^J- + 2 E E Qi t^^-' ^^-'l + 2 E l^ki^ "i^ + "i'] "^^ 



fc=2 



i=2 ;=i 



i=2 1=1 



2<k<l<n 



Using Lemma 5 i? is rewritten as 

n n n n n 

i? = 2E K,0fc] 9fe -2mE K,0fc] 5i - (i + m)E Kfc,0fe] +"^EE Kfc' [""'%]] 



fc=2 



fc=2 



fe=2 



fe=2 ;=2 
l^k 



-2E[0,«ifc]+2E[%"i^]^^- + 2EEQi"j' + 2 E [©feP^ifc + ^^n] 



Mfei 



fe=2 



J=2 



J=2 /=1 



2<k<l<n 



18 



and, hence, 

n n n n n 

R = -2m ^ [uifc, e^] 9 - (1 + to) ^ [w'lfe, e^] + TO ^ ^ [wife, [uK, e^J ] + 2 ^ [6, Wife] 

fc=2 fc=2 fc=2 1=2 fe=2 

n 

j=2 2<k<l<n 2<k<l<n 

Now, let us apply Lemma 4 which simplifies the above expression into the following 

n n n n n 

R = {l + m)Y,<k^k + (1 -"')E%^"ife - 2"^EE [Qfei'^ifc] "" - (1 + "^)E"ifeQfe 

fc=2 k=2 k=2 1=2 k=2 

l^k 



-(l + TO)^e^u'ife+TO^^[Mlfe,[uU,e^,-]] + Y^ [Q^,2{uiu+Uii)\ 

k=2 fe=2 1=2 2<k<l<n 

l^k 

n 

-2E©Xi. + 2 E i^k-QiXi 

j=2 2<k<l<n 



Ukl 



n n 



= -2to^^ [e^;,uife] ui; + TO^^ [uife, [ui/,ej,]] 

fc=2 /=2 fc=2 ;=2 

l^k l^k 

+2 E [e^p("u + "i/)]"fe^ + 2 E {^k-^iXi 



2<k<l<n 2<k<l<n 

From (|ll|l 9^, — 0; = [q — i^fc, 8^,;-] and, therefore, R is simplified to 



n n 



i? == "^EE ["l'="l''®fei] + 2 E [©M, (Wlfc + WiO"h] +2 ^ [i;i; -<;ife,6fc,] Uj 

fe=2 1=2 2<k<l<n 2<k<l<n 

l^k ^ ^ ~ ~ 

= 2 E [{^i~^k)u'f.i- {uik + uu)uki+muikUu,Qf,i\. 



2<k<l<n 

All derivatives with respect to xi of the term 
{q - Tfc) u'^.l - {uik + uii) Ukl + muikUii 

= TO (to + 1)^ I ( <; {Xi - Xl) -(, {Xi - Xk) ] <Pki - (pifc + P\i) Pki + pikpii 



are zero. Indeed, this follows from the fact that the first derivative is the addition theorem for the Weierstrass 
p- function. 

— I U{xi -xi) -i,{xi -Xk)j p'ki - (pife + pii)pki + Pikpii I 

= (plfc - pi/) p'ki - (p'lfe + p'li) pki + p'lkPll + Pikp'ii = 0. 
Therefore i? = 2 Y,2<k<i<n [('^i' ~ '^ifc) "w ~ (^ife + "lO "fe' + "^MlfcUl^, O^f] = and Lemma 6 is proven. 
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Proof of Theorem 1. 

We use Lemma 3 to reduce the commutator [/, H] to the expression l(T7)> and we show below that ((T7|l is zero. 
From ^ 



*=1 cre6(i;2t) 



therefore we can rewrite (|17f) as foUowing 



X — ^n 



Yx^j^didk,2uik 



.k=2 



*=1 o-6e(i;2t) 



El, hi+E E ^.e,Ja„25]u,, 



'j=2 



aee(i]:2t) 



1=1 

1^3 



( 



E %+ E E ^.0Jin^»^^'2^ 

2<i<i<n \ t=\ a^&(\if,2i) 



[Q'^]+E E [^.,-A.]0^+E E 

*=1 cree(i;2t) *=1 (Te6(i;2t) 



i=2 



X„ 



+E' 



J=2 



0i9,,2^ 



Uj7 



/=1 



j:'L,j: j: 



^j=2 



ELE j: 



^J=2 



cree(ii;2t) 



*-l <Te6(ii;2t) 
n 



x„a„2^ 



Ujl 



le<J 



0ja 



/=1 



Xn 



2<i<j<n 



Bi 



+ E E E [^.9.a„2i.,,]ejj^ 

2<i<]<n t=l ^ge(iii;2t) 



Xj9i,2^Ui; 



Z=2 



^Xjj,9i(9fc,2uii 



.fe=2 



From the formulae for the system with n = 2 we have 

[^{fc,i}'-A{fe,a] = -[dkdu2uki\, and [c^i^fc, 2uifc] = - [e{i,fc}, ff{i,fc}] 

and from Lemma 2 we have that for any c' G 6 (l; 2i) , i > 2, 

[X„.,-A,.] = -2 ^ [X,,\{k,i}dkduuki\ 



k<l 



and 



t=i ^Ge(<T'\{j};2t) 






Q-,a = 0. 
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Hence, [/, H] can be simplified to 



[I,H] = [Q,H]+J2 E [Qa,Ha]X,+ J2 E 

*=1 o-6e(i;2t) *=1 o-6e(i;2t) 



6^,2^. 



'ij 



jeer 



X„ 



+E 



J=2 



Qi^i'2^uj/ 






*-i creS(ii;2t) 



^j»dj,2 Y^ uji 



1=1 



x„ 



X£i9i,2^Ui; 



1=2 



E [0{i,fc}'-f^{i,fc}] ^ifc- 



fe=2 



This expression is equal zero. To verify this we should use relation p4|) . apply Lemma 6, which gives us 



and 



[e,i/] = -E 

J=2 



i=2 



0i^j'2E"j' 






^ [ej,,2(uife + Mij)]x 



{fc,i} 



2<fe<i<n 



e,j9„2^ 



Ujl 



1=1 



- E [Qakb'^i^ik + un)\X. 



{fe,i}, 



2<fe<i<n 



and equality 

[0{i,fc,O' ^{i,fc,o] = ~ [^1' 2 (uifc + ui/)] ^{fe,i} - [Q{i,i}dk, 2 (uifc + Ufc()] - [0{i,fc}5/, 2 (u^ + Ufe;)] , (22) 

The last equality can be obtained by direct calculation or using formulae for the three-particle case. Once all 
these substitutions are made it can be easily seen that all terms in l(T7| are canceled. This completes the proof of 
Theorem 1. 

7 Concluding remarks. 

A more general class of the deformed CM operators related to any Lie superalgebra has been recently introduced in 
|16j by Sergeev and Veselov, who found also recurrent formulae for the quantum integrals in trigonometric case for 
the classical series. In the elliptic case the integrability of these systems is an open problem. Within this approach 
the operator (^ considered in our paper corresponds to the Lie superalgebra sZ(n — 1, 1). The author hopes that 
the technique developed in this paper can be applied to more general deformed elliptic CM operators related to 
Lie superalgebra sl{n,p). 

In this relation it is worth mentioning that as one can see from our formulae ^ the parameters m = 1/Z, I — 
1,2,3... play a special role in the theory of the deformed CM operators. It is interesting to compare this with 
the results of the papers [IHI [T7j , where exactly these particular values of the parameter arise in relation with the 
strata in the discriminant variety consisting of polynomials having a root of multiplicity I. 
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